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Abstract 

We calculate the 1-loop renormalization of the fermion self-energy, all local 
fermion bilinears, as well as a set of extended bilinears which form a basis cor- 
responding to moments of the parton distribution functions. 

We use the overlap action for fermions and Symanzik improved action for gluons. 

Our results are presented as a function of the overlap parameter p and the 
parameters entering the Symanzik action. 



1 Introduction. 



In recent years, simulations in lattice QCD have employed an ever increasing number of 
"improved" actions for gluon and fermion fields. This practice is well motivated, since it 
leads to a reduction of 0(a) effects (a: lattice size), thus resulting in a faster approach 
to continuum physics. In the case of fermions, certain improved actions carry the crucial 
added advantage of preserving chiral symmetry on the lattice. 

The variety of lattice actions presently in use calls for the computation of a large num- 
ber of renormalization functions, before a comparison to physical results can be attempted; 
such functions often pertain to fermion composite operators, whose matrix elements are 
widely used to extract mass spectra, decay rates and transition amplitudes. 

A set of operators which are of particular interest are fermion bilinears; they can 
be classified according to the representations of the rotation group, and the global flavor 
group. Operators of this kind appear naturally in the operator product expansion (OPE), 
as applied to hadronic physics; in the Bjorken limit, involved in studies of deep inelastic 
scattering, the OPE for a product of hadronic currents takes the form: 

J(x)J(0) ~^C"' , ( a ; 2 )/ 1 ...^Oj^(0) (1) 

n,i 

The forward matrix elements of the local operators O appearing in this expansion 
are directly related to the moments of hadron structure functions. The dominant contri- 
bution in the expansion is given by operators whose twist (dimension minus spin) equals 
two, which in the flavor non-singlet case means the symmetric traceless operators pQ 

A a 

CW-Mn = ipj {fl D fil ...D fln} —i) (2) 

Of^^hi^D^-.D^^ (3) 

where A a are flavor matrices, and the curly brackets denote symmetrization over Lorentz 
indices. A proper renormalization is required for these operators, before one can relate 
their matrix elements, as extracted from numerical simulations, to physically observable 
quantities. 

In this work we employ Symanzik improved actions for gluons. We consider, in partic- 
ular, actions made up of closed Wilson loops with up to 6 links; they are the most widely 
used cases in numerical simulations. This class of actions is described by 3 independent 
parameters which are tuned either in perturbation theory (at tree level or beyond) or 
through simulation, so as to minimize 0(a) effects. Particular choices of action within 
this class, aside from the standard Wilson action, are: tree level Symanzik improved, 
"tadpole improved" Luscher-Weisz, Iwasaki and DBW2. 

For the fermion part of the action we employ Neuberger's overlap-Dirac operator. This 
action is characterized by one free parameter. The definition of this operator involves the 
inverse square root of X'X, where X is the Wilson-Dirac operator with a negative mass 
term; this operation is extremely costly in numerical simulations, and its viability depends 
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critically on the condition number of X*X. Thus, an improved action for gluons becomes 
all the more important, since it tends to decrease the condition number [2|. 

The quantities whose renormalizations we set out to compute, to 1-loop order in lattice 
perturbation theory, are: 



a) The quark propagator 

b) Local fermion bilinears. 



(Mx)Mv)) (4) 



O x = i>(x)T x ij(x) (5) 

c) Operators which measure the first moment of quark momentum distributions. 

V2 = $"/ {1 D 4} Tp (6) 

0' V2 =i> lA D^-~J2i>liD^ (7) 
6 i=i 

d) Operators which measure the first moment of quark helicity distributions. 

Oa 2 = i)l{ilhD A} il) (8) 

1 3 _ 

0'a 2 = ^7475^4^ - o ^HlhDiil) (9) 
6 i=i 

The symmetrized, covariant derivatives D are defined as: D = D — D , with: 

D fl ip(x) = — U(x, /i) ip{x + afx) — U\x — ajl, /i) ip{x — afi) (10) 

and similarly for the left derivative D . 

Our calculations provide a cross check and an extension of results appearing in Refs. 
IHH IH E]. Work is in progress using operators measuring the second moment of quark 
momentum and helicity distributions. 

We have performed our calculations for a rather wide selection of values for the pa- 
rameters entering the gluon and fermion actions, corresponding to present and future 
simulations. In particular, our results scan the whole range of admissible values for the 
overlap parameter. Results for different choices of parameter values can be provided by 
the authors upon request. 

A synopsis of some of our results can be found in Ref. jSj. 
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2 Calculational setup. 



We denote the lattice action by 



S — Sg + S_ 



'ID 



where Sq is the gluon action, and Sp is the fermion action. The gluon action we consider 
is written in standard notation: 



Sg = — \(c Q Trllpi + C! ^ TrU rtg 

9 plaquette rectangle 

+C 2 TrU chr + C 3 TrU plg) 

chair parallelogram 



(12) 



where U v \ is the standard plaquette, while the remaining U's cover all possible closed 
loops containing up to six links, as indicated in Fig. 1. 



plaquette 



rectangle 



chair = 




Figure 1: Closed loops corresponding to the four terms in the gluon action (|12|) . 

The coefficients cq, c±, c 2 , C3 satisfy the normalization condition: 

c + 8ci + 16c 2 + 8c 3 = 1 (13) 

As in Ref.|4j, we have used the values of Co, ci, C2, C3, shown in Table 1, where "Plaquette" 
is the standard Wilson action for gluons, "Symanzik" is the tree-level improved action 
[7j and TILW is the tadpole improved Luescher- Weisz action jSJ El . Two further choices 
which we have considered for q are the Iwasaki ^0] and DBW2 [llj actions. It is worth 
noting that, at the 1-loop level, these coefficients enter the calculation only through the 
propagator and not through any vertices; given that the propagator does not depend 
separately on c 2 and C3, but only on their sum, we have considered only parameter values 
with c 2 = 0. 

The action for massless overlap fermions is given by ^2] 

s^EEfi^nW ( 14 ) 

/ x >y 

with 

D N = p[l+X(X t X)-5] (15) 
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Act ion 




Co 


Cl 


C2 


C3 


Plaquette 




1 n 
l.U 


n n 
U.U 


n n 
U.U 


n n 

U.U 


Symanzik 




1.66666 


-U.Uooooo 


U.U 


U.U 


TTT,W1 13 — 8 

J. -LU V V X LJ (J 


fin 


9 31 681 


-0 1 ^1 791 





-001 98098 


TILW2 (3 = 8 


.45 


2.34602 


-0.154846 


0.0 


-0.0134070 


TILW3 (3 = 8 


.30 


2.38698 


-0.159128 


0.0 


-0.0142442 


TILW4 (3 = 8 


.20 


2.41278 


-0.161827 


0.0 


-0.0147710 


TILW5 (3 = 8 


.10 


2.44654 


-0.165353 


0.0 


-0.0154645 


TILW6 (3 = 8 


.00 


2.48917 


-0.169805 


0.0 


-0.0163414 


Iwasaki 




3.648 


-0.331 


0.0 


0.0 


DBW2 




12.2688 


-1.4086 


0.0 


0.0 



Table 1: The values of coefficients cq, c±, C2, C3 . 

and: X = D w — p. Here .D,,, is the massless Wilson-Dirac operator with r = 1, and p is 
a free parameter whose value must be in the range < p < 2 in order to guarantee the 
correct pole structure of D N . 



3 Results. 

Quark propagator 

Let us consider the massless quark propagator SV first. In momentum space, the 
inverse of can be written as: 

with C F = (N 2 — 1)/2N, and £i(a,p) = log(a 2 p 2 ) + 6 S (Feynman gauge). In Eq. (fTo^) . 
the factor in parentheses is the quark wave function renormalization; in the MOM scheme 
we have: 

ZF OM (a,fj) = 1 - f^f£i(^) (17) 

The diagrams that contribute to one loop order are shown in Fig. 2. Our results for 
are shown in Fig. 3 and, on a different scale, in Fig. 4. One observes a divergence of 
this quantity for p — > 0, which can be expected, since the fermion propagator behaves 
singularly in this limit. Numerical values of 6s are listed in Table 2 for typical values of 
p: p = 0.6, 1.0, 1.4. 
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Figure 2: One-loop diagrams contributing to the quark self-energy. 
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Figure 3: 6g as a function of the overlap parameter p, for different gluon actions. 
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Figure 4: 6s as a function of the overlap parameter p, shown on a different scale. 



Quark bilinears 

Let us consider local operators of the form 



Ox = i){x)T*i){x) (18) 

where X = S, P, V, A, T, and V s = 1, T p = 7s , T v = 7m , T a = 7^75, T T = . We 
denote the amputated Green's function of the operator Ox by A x . The final results are 
(in Feynman gauge): 



A s > p = {1 )75 } + t^T [-4MaY) + b S ,p]{hl5} 



9 2 C 



16tt 



A 



V,A 



9 2 C 



{7^, 7^75} + T^r[ _ 7/,(^^(a p ) + b ViA ) + 2— ]{1, 75} 



IQtt 



Mm. 



(19) 

(20) 
(21) 



Where 65 = 6p and 6y = 6a- Among the diagrams appearing in Fig. 5, only the first one 
(Vertex) contributes to A x . 
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Figure 5: One-loop lattice Feynman diagrams contributing to the amputated Green's functions 
of Ox and Onu,i...u m - From left to right: vertex, operator tadpole, left sail and right sail. 

Our results provide a cross check and an extension of results appearing in Refs. [31 H] . 
They are shown in Figures 3+4 (6s), 6 (6s,p), 7 (bv,A), 8 (br) for typical values of the 
overlap parameter p. The same quantities are also listed in Tables 2, 3, 4, 5, respectively, 
for the representative choices: p = 0.6, 1.0, 1.4. In all cases which can be compared, our 
results agree with those of Refs. [3J Hj. There is one exception, regarding the values of 
for cases other than the standard plaquette \?^ aq : In these cases, our results for b T — bj, aq 
have the opposite sign compared to Ref. j3]. 
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Figure 6: bs t p as a function of the overlap parameter p, for different gluon actions. 
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Figure 7: by^A as a function of the overlap parameter p, for different gluon actions. 
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Figure 8: &t as a function of the overlap parameter p, for different gluon actions. 



8 



First moment of quark distributions. 



We have calculated the 1-loop renormalization coefficients of the operators Ov 2 ,0'y 2 
which are the symmetric off- diagonal, and the traceless diagonal parts, respectively, of 
the operator: 

= in^Dvyb (22) 

The axial counterparts Oa 2 ,0' a of the above operators renormalize in the same way. 
Using the notation of Ref . [5] , we find for the amputated Green's function of 0^ v in the 
Feynman gauge: 



9 2 C, 
16tt 2 



r,5. 



{-log(a p)+b) 1{^Pu} + b' 5^ v p v - 



4 y^Vu 

3 9 J 



(23) 



The quantities b, b' correspond to (pi + 62), &4 of Ref. respectively. All diagrams of 
Fig. 5 contribute in this case. The rational coefficients 5/3 and -4/3 in Eq. (J2H|) check 
with those of Ref. [3]. The values of b and b' are shown in Figures 9 and 10; they are in 
agreement with the quantities bi + b 2 and 64, respectively, as defined in Ref. [3], for the 
choice of parameters considered in that reference. In Tables 6 and 7 we list our results 
for b and b', respectively, with the choice of values for p: p = 0.6, 1.0, 1.4. 
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Figure 9: b as a function of the overlap parameter p, for different gluon actions. 
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Figure 10: b' as a function of the overlap parameter p, for different gluon actions. 

All 1-loop calculations presented in this work were performed using our Mathemat- 
ica package for lattice perturbation theory. Numerically inverting the gluon propagator 
turned out to be less CPU consuming, as compared to using an analytic expression for 
the inverse. We used a rather wide spectrum of values for the Symanzik parameters q ; 
given that the dependence of 1-loop results on c 2 and c 3 is only through their sum c 2 + c 3 , 
we have chosen c 2 = without loss of generality. For the overlap parameter p we have 
explored in detail its full spectrum of admissible values. Results for further values of the 
parameters can be provided by the authors upon request. 

We observe a rather smooth dependence of our results on p, and this fact allows for a 
safe interpolation to any additional p-values in the allowed range < p < 2. A possible 
exception is the behaviour of &x near the endpoint p — * 0, which can be expected given 
the incorrect pole structure of the fermion operator in that limit. The dependence on p 
is even less pronounced for the DBW2 action. 

Results for the various choices of parameters for the Tadpole Improved Liischer-Weisz 
(TILW) action are in many cases almost indistinguishable from one another, due to the 
proximity in these parameter values; consequently, it is straightforward to interpolate 
accurately our results to other sets of values which were not considered in this work. 

It is worth noticing that both b and b' are noticeably suppressed in the case of the 
DBW2 action. A similar effect has been observed in the perturbative renormalization 
of both the topological charge and susceptibility, using this action 13J. This fact is an 
indication in favour of this action for purposes of numerical simulation. 

Some related quantities of interest for the study of structure functions are the second 
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moments of quark momentum and helicity distributions, as well as non-forward matrix 
elements of various fermion bilinears. We expect to address these issues in a future report. 
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TABLES 



A c*^ l r\ n 






n — 1 n 


n — 1 4 


JT ldL[U.t:L It; 








93 7fi^Q7Q1 Q^lVzL 


o y iiidiizjij^. 






3D ^D4Q1 Q93QQf3^ 


1 Q 731 QQ c i9Q4 c il'9 


TILWl /3 = 8 


60 


-48.315432542(1) 


-26.5273492659(2) 


-17.5090552592(4 


TILW2 = 8 


45 


-48.020432621(1) 


-26.3832716142(5) 


-17.4287897913(1 


TILW3 = 8 


30 


-47.614837428(1) 


-26.1852273116(3) 


-17.3184864381(3 


TILW4 /3 = 8 


20 


-47.363843382(1) 


-26.0626992727(4) 


-17.2502579516(4 


TILW5 (3 = 8 


10 


-47.040843009(1) 


-25.9050530413(4) 


-17.1624938296(4 


TILW6 = 8 


00 


-46.641145138(1) 


-25.7100225495(5) 


-17.0539447038(2 


Iwasaki 




-37.911348699(2) 


-21.454735885(1) 


-14.680753624(5) 


DBW2 




-19.673287561(6) 


-12.675353551(3) 


-9.840831462(1) 



Table 2: The contribution 6 S to the self energy for various actions. 



Action 




p = 0.6 


p = 1.0 


p = 1.4 


Plaquette 




8.9798281638(5) 


10.26190085663(6) 


11.34335740638(2) 


Symanzik 




8.6712333358(5) 


9.87030941138(2) 


10.85599048485(3) 


TILWl (3 = 8 


.60 


8.4347849739(6) 


9.57373721242(1) 


10.4927473931(1) 


TILW2 (3 = 8 


.45 


8.4249104780(5) 


9.56141189336(1) 


10.4777493701(1) 


TILW3 (3 = 8 


.30 


8.4111618345(5) 


9.54425880565(2) 


10.45688939118(1) 


TILW4 (3 = 8.20 


8.4025516764(5) 


9.53352128912(2) 


10.44383891932(1) 


TILW5 (3 = 8 


.10 


8.3913573861(6) 


9.51956660551(1) 


10.4268869200(1) 


TILW6 (3 = 8 


.00 


8.3773207579(2) 


9.50207732589(2) 


10.4056548274(1) 


Iwasaki 




7.989846185(1) 


9.02297610282(2) 


9.82968352735(1) 


DBW2 




6.346306461(1) 


7.0639464328(2) 


7.571097190(2) 



Table 3: The contribution bs,p to A s,p for various actions. 



Action 




p = 0.6 


p = 1.0 


p = 1.4 


Plaquette 




6.3204575633(1) 


6.33125072110(1) 


6.34801586915(1) 


Symanzik 




6.3141428350(1) 


6.32252999501(1) 


6.3354222684(1) 


TILWl (3 = 8 


.60 


6.3107172748(1) 


6.31780175061(1) 


6.32862040074(3) 


TILW2 (3 = 8 


.45 


6.3105909672(1) 


6.31762724546(3) 


6.32836942659(2) 


TILW3 (3 = 8 


.30 


6.3104170577(1) 


6.31738695411(2) 


6.32802384701(4) 


TILW4 (3 = 8 


.20 


6.3103092794(1) 


6.31723802434(3) 


6.32780966499(3) 


TILW5 (3 = 8 


.10 


6.3101704948(1) 


6.31704624086(2) 


6.32753386701(3) 


TILW6 (3 = 8 


.00 


6.3099985210(1) 


6.31680857518(1) 


6.32719209974(1) 


Iwasaki 




6.3057833598(1) 


6.31095730979(1) 


6.31875098630(4) 


DBW2 




6.2975701799(1) 


6.2996002928(1) 


6.30255799548(6) 



Table 4: The contribution b v ,A to A V ' A for various actions. 
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Action 




p = 0.6 




p = 1.0 


p = 1.4 


Plaquette 




4.7673340248(1) 


4 


35436733681(1) 


4.0162353498(7) 


Symanzik 




4.86177933052(5) 


4 


4732701854(1) 


4.16189952383(6 


TILWl (3 = 8 


.60 


4.9360280383(1) 


4 


56582325908(4) 


4.27391139850(6 


TILW2 = 8 


.45 


4.93915112696(5) 


4 


56969902605(5) 


4.27857610740(2 


TILW3 = 8 


.30 


4.94350212841(2) 


4 


57509633367(4) 


4.28506866099(4 


TILW4 /3 = 8 


.20 


4.94622847715(4) 


4 


57847693279(6) 


4.28913324224(3 


TILW5 /3 = 8 


.10 


4.9497748610(1) 


4 


58287278284(4) 


4.29441617825(3 


TILW6 p = 8 


.00 


4.95422443865(1) 


4 


58838565488(6) 


4.30103785269(2 


Iwasaki 




5.0777624156(2) 


4 


7402843779(1) 


4.48177346976(3 


DBW2 




5.6146580874(6) 


5 


378151584(1) 


5.213044934(1) 



Table 5: The contribution br to A T for various actions. 



Action 




p = 0.6 


p = 1.0 


p= 1.4 


Plaquette 




-15.7863271766(1) 


-15.6250817532(1) 


-15.4858619898(1) 


Symanzik 




-11.9937499031(2) 


-11.8664388808(2) 


-11.7640746811(2) 


TILWl p = 8 


.60 


-9.7160722504(2) 


-9.6095867774(2) 


-9.5287008568(2) 


TILW2 (3 = 8 


.45 


-9.6309328038(2) 


-9.5252267216(1) 


-9.4451255929(1) 


TILW3 p = 8 


.30 


-9.5135691580(2) 


-9.4089372737(1) 


-9.3299153795(1) 


TILW4 (3 = 8 


.20 


-9.4407612246(2) 


-9.3367955747(1) 


-9.2584416246(1) 


TILW5 p = 8 


.10 


-9.3468642074(4) 


-9.2437578051(4) 


-9.1662637020(4) 


TILW6 p = 8 


.00 


-9.2303490445(1) 


-9.1283086414(1) 


-9.0518787415(1) 


Iwasaki 




-6.5572678902(2) 


-6.4778880567(1) 


-6.4232808744(1) 


DBW2 




0.185440914(4) 


0.211859722(4) 


0.222852646(3) 



Table 6: The contribution b to A,,„ for various actions. 



Action 




p = 0.6 


p = 1.0 


p= 1.4 


Plaquette 




0.3497088259(1) 


0.4340050439(1) 


0.5267146951(2) 


Symanzik 




0.2719708888(4) 


0.3367316099(2) 


0.4080895543(1) 


TILWl p = 8 


.60 


0.2307195514(3) 


0.2858316439(2) 


0.3469216320(2) 


TILW2 (3 = 8 


.45 


0.2290786112(2) 


0.2838197180(1) 


0.3445199656(2) 


TILW3 p = 8 


.30 


0.2268014376(4) 


0.2810293893(1) 


0.3411911844(3) 


TILW4 p = 8 


.20 


0.2253789601(2) 


0.2792873534(1) 


0.3391142249(3) 


TILW5 p = 8 


.10 


0.2235391956(1) 


0.2770353882(1) 


0.3364306054(3) 


TILW6 p = 8 


.00 


0.2212413920(1) 


0.2742245511(1) 


0.3330831698(2) 


Iwasaki 




0.1277382408(3) 


0.1614096997(1) 


0.2010432525(2) 


DBW2 




-0.116651374(1) 


-0.118539399(1) 


-0.111845323(1) 



Table 7: The contribution b' to A,,„ for various actions. 
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